Tunneling current spectroscopy of a nanostructure junction involving multiple energy 
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A multi-level Anderson model is employed to simulate the system of a nanostructure tunnel 
junction with any number of one-particle energy levels. The tunneling current, including both shell- 
tunneling and shell-filling cases, is theoretically investigated via the nonequilibrium Green's function 
method. We obtain a closed form for the spectral function, which is used to analyze the complicated 
tunneling current spectra of a quantum dot or molecule embedded in a double-barrier junction. We 
also show that negative differential conductance can be observed in a quantum dot tunnel junction 
when the Coulomb interactions with neighboring quantum dots are taken into account. 
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In the last decade, the tunnehng current spectra of 
nanostructure junctions such as single-electron transis- 
tors (SETs)-^"^ and molecule transistors (MTs)^~^ have 
been extensively studied due to their important appUca- 
tions in quantum computing, quantum communication, 
and ultrahigh-density IC circuits. Although the single- 
particle energy levels of such systems can be obtained 
by ah inito or semiempirical methods^" it is still dif- 
ficult to model the full tunneling current spectra due to 
the presence of Coulomb blockade and the uncertainty in 
dot shape, size and position. 

The characteristics of tunneling current of SETs or 
MTs can be classified into either the shell-tunneling 
case (with no charge accumulation) or the shell-filling 
case (with charge accumulation). In the shell-tunneling 
case, the electron Coulomb interactions are suppressed. 
Consequently, current spectra directly reveals the one- 
particle resonance energy levels. In the shell-filling case, 
the tunneling current spectra become much more com- 
plicated due to the presence of electron Coulomb inter- 
actions. Consequently, the charging energies and energy 
levels can not be readily and directly obtained from ex- 
perimental data. Both physical parameters are crucial 
not only in the optimization of tunneling devices, but also 
for the understanding of fundamental physics of minia- 
turized nanostructures. Thus, it is desirable to derive a 
simple formula that allows one to determine the multiple 
energy levels involved and the associated Coulomb inter- 
actions from the measured tunneling current spectra. 

In this letter, we derive a simple yet general analytic 
expression that is applicable for modeling the tunnel- 
ing current through any nanostructure (including QDs 
and molecules) involving multiple energy levels (i.e. the 
ground state plus all the excited states of interest). This 
is done by solving a multi-level Anderson model via non- 
equilibrium Green's function technique, which has been 
extensively used for investigating the Coulomb block- 
ade and Kondo effect on the tunneling current through 
the ground state of a single QD.^^ We find that when 
more than one level in a QD/molecule are involved in 
the charge transport, the average two-particle occupation 
numbers play an important role in determining whether 



a given energy channel is close or open to the tunneling 
current. 

The system of an isolated nanostructure embedded in a 
double-barrier junction can be described by the following 
Hamiltonian i3-i6. 

+ Y ^k,(l,e4,a'^k,a,/3, (1) 

where a\ ^ ^ {ak.a-./s) creates (destroys) an electron of mo- 
mentum k and spin a with energy ek in the P metallic 
electrode. (i| ^ (rf^.o-) creates (destroys) an electron inside 
the nanostructure with orbital energy Ei. Uij describes 
the intra- or inter-level charging energies, and Vk,i3,i de- 
scribes the coupling between the band states in the con- 
tacts and the confined states in the nanostructure. 

Using Keldysh Green's function technique"'^"'^'^^, we ob- 
tain the tunneling current through a nanostructure. 
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where /l = /(e - ^l) and Jr = f{e - (ir) are the 
Fermi distribution functions for the left and right elec- 
trode, respectively. G^^(e) is the retarded Green's func- 
tion. The chemical potential difference between these 
two electrodes is equal to the applied bias eVa- Ti^l and 
^i,R [^i,f3 = I]k \ Ve.p,k\'^S{e - ek)] denote the tunneling 
rates from the nanostructure to the left (source) and right 
(drain) electrodes, respectively. For simplicity, these tun- 
neling rates will be assumed to be energy independent. 
Therefore, the calculation of tunneling current is entirely 
determined by the spectral function, A{e) — ImG\ ^{e). 

Rigorous solution to ^ in the Coulomb blockade 
regime can be obtained by solving a hierarchy of equa- 
tions of motion, which relate G\ ^ to two-particle Green's 
functions, which can be solved due to the termination of 
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the series at some finite order. After tedious algebra, we 
obtain a neat closed form (by keeping the two-particle 
correlation functions associate with the same level) 
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where n denotes the number of one-particle energy levels 
considered for the nanostructure. !!„ denotes the sum 
of Coulomb interactions seen by a particle in level £ due 
to other particles in configuration m, in which each level 
j{j 7^ ^) can be occupied by zero, one or two particles, pm 
denotes the probability of finding the system in configu- 
ration TO. For a two-level (n — 2) system {£ ^ j), we have 



three configurations with pi 



l-(7V,-<, + iV,-_.) + 



(the probability with no particle in level j), 
P2 — V = Nj^rj + ^j,-cr — '^i'rij.a-nj^-cr) (the probability 
with one particle in level j), and p^ = = (rij^a-nj^-a) 
(the probability with two particles in level j). Mean- 
while, El = 0, 112 = Uej and lis — 2Uej. Ugj denotes 
the interlevel Coulomb interaction between an electron 
at level £ and the other at level j. For a three-level 
system {£ ^ j ^ /), there are nine (3 x 3) configura- 
tions, and the probability factors become pi — a^a^ , 
P2 — V , P3 — a^V , Pi — c^a^ , p^ — , pe — , 
Pt = d'V , ps = c' V , and j3g ^ d> cJ . Interlevel Coulomb 
interaction factors are Hi = 0, Ila = Uij, Ha = U^ji, 
n4 = 2Ui,j, Hs = 2Ui,j,, He = Ugj+Uij-, Ur = 2Uij + Uij,, 



Us ~ 2Uej' + Uej, and Ilg = 2Uij + 2Uiji. Based on 
these simple rules , the probability factors and interlevel 
Coulomb interaction factors for any number of energy 
levels can be similarly determined. We see that the prob- 
ability of finding the system in each configuration is de- 
termined not only by the average one-particle occupation 
numbers but also by the average two-particle occupation 
numbers. N^^„ in Eq. (3) can be obtained by solving the 
following equations self-consistently 
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Both iVf.cr and Ni^i = {nf^a-ne,-a-) sue limited to the range 
between and 1. 

To illustrate the usefulness of our theory, we calcu- 
late the tunneling current spectra of an isolated CdSe 
QD sandwiched between an STM tip (left lead) and 
a conducting substrate (right lead) as studied experi- 
mentally in Ref. We consider the case of three en- 
ergy levels. Because the tip-substrate is biased at T4, 



the bare energy levels of Ei in the dot are changed to 
El + aeVa . The factor a = Ll/{Ll + Lr) is determined 
by the distances from the QD center to the tip [Ll) 
and the substrate (ii?), respectively. Here a — 0.61, 
which is determined from the separation between the 
first peak at negative bias and that at positive biasJ^ 
Other physical parameters can also be determined by 
comparing Eq. (3) with the peak positions of the tun- 
neling spectra observed in Ref. We obtain the fol- 
lowing. The chemical potentials of both electrodes are 
0.78 eV below the ground state level Ei at zero bias. 
E2-Ei= 0.236 eV, E3-Ei = 0.456 eV, Ui = 0.137 eV, 
Ui2 = U21 = 0.122 eV, U2 = 0.07 eV, U3 = 0.06 eV, 
Ui3 = U31 = 0.1 eV, and U23 = U32 = 0.04 eV. These 
intra- and inter-level Coulomb energies are reasonable for 
the 3 nm diameter CdSe QD considered here. We note 
that the p-like level in a spherical QD is sixfold degen- 
erate (including spin). However, the coupling strength 
between the tip (or substrate) and the p^- or pj^-like or- 
bital is weak. Therefore, only px-like orbital has been 
included, which is labeled by E2, while E3 denotes one of 
the d-like orbitals that is strongly coupled to the leads. 

Occupation numbers Ni and Ni^i are obtained by solv- 
ing the coupled equations Eqs. (4)- (6). Once they are de- 
termined, we can calculate the tunneling current by sub- 
stituting Eq. (3) into Eq. (2). Figure 1 shows the differ- 
ential conductance, dJ/dVa as a function of applied bias 
for various ratios of Tl/Th. Curve (c) exhibits the max- 
imum number of resonant channels, which correspond to 
the poles of the retarded Green's function in Eq. (3). The 
first eight peaks of curve (c) (with energy less than £^3) are 
ei = El, €2 = Ei+Ui, £3 = E2, e4 = E1 + U1 + U12, £5 = 
E2 + U2, ee ^ E2 + C/12, €7 = El + Ui + 2Ui2, and 
es ^ E2 + U2 + Ui2- The charging energy of the ground 
state, Ui can be determined by the difference in bias be- 
tween the first two peaks. The two channels, e = E1 + U12 
and e = Ei+2Ui2 are prohibited because those resonance 
energies are below E2 (hence N2 — and iV22 = when 
the applied bias aligns with these channels). The limit 
^l/^r ^ 1 corresponds to the so called " shell-tunnefing 
case".^^ According to Eqs. (4) and (5), and Ni^i are 
determined by the factor V = T£^L/(r£^L + ^e,R)- There- 
fore, Ni and Ni^i approaches zero in this limit. Con- 
sequently, the peaks associated with the resonance lev- 
els involving the intralevel and interlevel Coulomb inter- 
actions will be suppressed, since the carrier will tunnel 
out much faster than it can be fed. The peaks associ- 
ated with the single-particle levels, Ei, E2 and E3 will 
be the only resonance channels allowed in the limit of 
ri/Ffl <JC 1 as seen in curve (e). In the opposite limit, 
r L I^R 3> 1 we have the so called " shell- filling case" [see 
curve(a)]. In this case, both iVi and A^n approach 1 
when Va > Ei + Ui. Consequently, resonance channels 
related at E2 and E2 + U2 are suppressed while the res- 
onances at E2 + 2[/i2 and E2 + U2 + 2Ui2 are enhanced. 
This indicates that the bare energy levels are renormal- 
ized if charges reside in the QD. Similarly for E3 related 
channels. 
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Next, we compare our theoretical prediction for the 
tunnehng current spectra with those measured in Ref . [13, 
where the observed differential conductance peaks are 
broadened. The broadening is mainly due to coupling 
to nearby QDs since the experiment was performed at 
T = 4K. To take into account the above broadening 
effect, we replace each Lorentzian function appearing in 
the differential conductance by a Gaussian function of 
the form exp{- ^^^^^j-' }/{p,V2t:). f,, and pi de- 
note the peak strength, resonance energy, and broaden- 
ing width, respectively. Figure 2 shows the predicted 
differential conductance as a function of applied bias at 
zero temperature. We have used the following parame- 
ters Tl.i = 1 meV, Tr^i = 0.15 meV, r^^a = 3 meV, 
rfl,2 = 0.6 meV.rL,3 = i.5 meV, Tr.s = 0.375 meV. The 
variation in tunneling rates for different levels reflects the 
difference in wave functions and the bias-dependent bar- 
rier height. For simplicity, we assume pi — p + Ti, where 
p = 35 meV characterizes the broadening due to coupling 
to nearby QDs. Both the positions and relative strengths 
of these peaks are in very good agreement with the exper- 
imental measurement reported in Ref. il7i . If we remove 
the extrinsic broadening by setting p = 0, the differential 
conductance spectrum exhibits more resonance channels 
as shown in the lower part of the figure, which can not 
be resolved experimentally. 

To illustrate the significance of the two-particle occu- 
pation number, we apply our theory to the case of a cou- 
pled dot embedded between two leads, where two QDs 
(dot A and dot B) are coupled by interdot Coulomb in- 
teractions, even though the inter-dot tunneling is weak. 
It is assumed that the right lead is closer to dot A than to 
dot B. Thus, the tunneling rate Tr b is ten times smaller 
than Tfi^A, and the total current is dominated by the cur- 
rent through dot A {Ja)- Figure 3 shows the total tunnel- 
ing current, J = Ja + Jb for various strengths of the in- 
terdot Coulomb interaction (which may, for example, cor- 
respond to different separations between the two dots). It 
is seen that negative differential conductance (NDC) can 
occur due to the interdot Coulomb interaction. The NDC 
characteristic of tunneling current can be understood 
mainly from the feature of J a, since J a ^ Jb- NDC 
first occurs at Eb as dot B becomes filled, which leads to 
a suppression of the current through the channel Ea, and 
the current resumes as the channel Ea + Uab opens up 
when the bias further increases. For comparison, we also 
plot in Fig. 3 the total current, J a + Jb obtained with 
the approximation, < nj^a-nj,-a- >— Nj^^Nj^^^ (dash- 
dotted curve) . The difference between this curve and the 



solid curve shows the significance of treating the two- 
particle occupation number correctly. In particular, we 
see in the dash-dotted curve ( Jq + Jb) that the resonance 
channel at e = Ea + "2X1 ab is allowed even though the 
applied bias is below Eb + Ub (i-e. dot B remains singly 
charged). Such a nonphysical behavior in Jq -|- Jb arises 
due to the approximation Nbb ~ NbNb 7^ used, while 
if treated correctly Nbb should vanish under this condi- 
tion. We also notice that the tunneling current increases 
only slightly when applied bias overcomes the charging 
energy of dot A (marked by Ea + Ua) since we are con- 
sidering a shell- filling case for dot B {Tl,b/^r,b = 10) 
here {Nb ~ !)■ Increasing the applied bias beyond the 
resonance channel at Ea + Ua + Uab leads to a significant 
jump in the tunneling current. However, this tunneling 
current is suppressed again at Eb + Ub + Uab because 
dot B is now filled with two charges. 

To further analyze the NDC behavior due to interdot 
Coulomb interaction, we plot in Figure 4 the differen- 
tial conductances for various ratios of ^l,b/^r,b- Curve 
(a) corresponds to the derivative of the curve shown in 
Fig. 3. The NDC behavior is quite apparent. However, 
the NDC behavior disappears in curves (b) and (c) when 
dot B is not under the shell-filling condition, even though 
the interdot Coulomb energy remains the same. When 
dot B is in the shell tunneling limit (F^.s ^ ^r.b), J a 
becomes unaffected by dot B, and the tunneling current 
[curve (c)] exhibits mainly two resonance channels at Ea 
and Ea + Ua plus a weaker peak at Eb + Uabi which 
is contributed from Jb- As illustrated in Figs. 3 and 
4, when one employs an SET as a charge detector, the 
interdot Coulomb interactions due to neighboring QDs 
and the tunneling rate ratios, Fl/F/j should be carefully 
considered in the design. 

In conclusion, a closed form expression for the spec- 
tral function for nanostructure tunnel junctions involv- 
ing multiple energy levels has been derived and it has be 
employed to analyze the complicated spectra of the tun- 
neling current through a realistic single QD. Our present 
theory also provides a useful guideline for the design of 
SET with coupled quantum dots. We show that when the 
interdot hoping effect is weak but the interdot Coulomb 
interactions are strong, NDC behavior can occur and it 
depends sensitively on the tunneling rate ratios for in- 
coming and outgoing electrons in each QD. 
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Figure captions 

Fig. 1 Differential conductance as a function of applied 
bias at for various ratios of Ti,/Tr: (a) 10, (b) 5, (c) 1, 
(d) 0.4, and (e) 0.2. Tr is fixed at 1 meV. 

Fig. 2 Differential conductance as a function of applied 
bias with and without the extrinsic broadening effect due 
to coupling to nearby QDs. For better display, the broad- 
ened spectra has been multiplied by a factor 5. 

Fig. 3 Tunneling current {Ja+Jb) as a function of bias 
for various strengths of interdot Coulomb interaction, 
Uab- The parameters used are: Ea = Ep + 108 meV, 
Eb = Ea + 10 meV, Ua = HI meV, Ub = 108 meV, 
Tl,a = Tr,a = Tl^b = 1 meV, and Tr^b = 0.1 meV, 
and a = 0.61. Current is in units of Jo=e*(meV)/h. 
Dot-dashed curve is for Ja + Jb with Uab = 30 meV. 

Fig. 4 Differential conductance as a function of ap- 
plied bias for various tunneling rate ratios o{Tl,b/^r,b- 
Uab = 30 meV. (a) Tl.b = 1 mcV and Tr_b = 0.1 meV. 
(b) rL,B = 1 meV and Tr^b = 1 meV. (c) F^^b = 
0.1 meV and Tr^b = 1 meV. 
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